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Abstract 

The scenario of gravitino dark matter with broken i?-parity naturally reconciles 
three paradigms that, albeit very well motivated separately, seem to be in mutual 
conflict: supersymmetric dark matter, thermal leptogenesis and standard Big 
Bang nucleosynthesis. Interestingly enough, the products of the gravitino decay 
could be observed, opening the possibility of indirect detection of gravitino dark 
matter. In this paper, we compute the positron and the antiproton fluxes from 
gravitino decay. We find that a gravitino with a mass of ~ 150 GeV and 
a lifetime of r3/2 ~ 10^^ s could simultaneously explain the EGRET anomaly 
in the extragalactic diffuse gamma ray background and the HEAT excess in the 
positron fraction. However, the predicted antiproton fiux tends to be too large, 
although the prediction suffers from large uncertainties and might be compatible 
with present observations for certain choices of propagation parameters. 
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1 Introduction 



Models with local supersymmetry predict the existence of a particle with extremely 
weak interactions: the gravitino. In contrast to the supersymmetric partners of the 
Standard Model particles, whose masses are expected to lie in the electroweak domain, 
the gravitino can have a mass ranging between a few eV and several TeV without 
conflicting with any laboratory experiment. Therefore, the gravitino can very naturally 
be the lightest supersymmetric particle (LSP), and if it is sufficiently long-lived, it could 
constitute the dark matter of the Universe [1]. 

Gravitinos were produced in the early Universe by scatterings in the thermal plasma, 
but did not subsequently annihilate due to their extremely weak interactions. There- 
fore, a relic population of gravitinos is expected in the present Universe with a density 
given by [2] 

where Tr is the reheating temperature of the Universe, rriz/i is the gravitino mass and 
m~ is the gluino mass. In predicting the relic abundance of gravitinos, the main uncer- 
tainty arises from our ignorance of the thermal history of the Universe before Big Bang 
nucleosynthesis (BBN) and in particular of the reheating temperature after inflation. 
However, we have strong indications that the Universe was very hot after inflation. 
Namely, the discovery of neutrino masses about ten years ago provided strong support 
to leptogenesis as the explanation for the observed baryon asymmetry of the Universe. 
This mechanism can reproduce the observed baryon asymmetry very naturally if the 
reheating temperature of the Universe was above 10^ GeV [3]. Therefore, following 
Eq. (1), the gravitino could constitute the dark matter if 777.3/2 ^ 10 GeV for a gluino 
mass m~ ~ 1 TeV, which is consistent with the assumption that the gravitino is the 
lightest supersymmetric particle. 

Remarkably, the conjectures of a reheating temperature of the Universe larger than 
10^ GeV and a gravitino mass larger than a few GeV can naturally solve two of the most 
long-standing problems in cosmology: the nature of the dark matter and the origin of 
the baryon asymmetry of the Universe. Nevertheless, this picture is not exempt from 
problems. If i?-parity is exactly conserved, the next-to-LSP (NSLP) can only decay 
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gravitationally into gravitinos and Standard Model particles with a lifetime 




(2) 



Then, the NLSP is typically present during and after Big Bang nucleosynthesis, jeop- 
ardizing the successful predictions of the standard nucleosynthesis scenario. This is in 
fact the case for the most hkely candidates for the NLSP: the lightest neutrahno and 
the right-handed stau (or more generically, any negatively charged particle, such as the 
chargino). More precisely, when the NLSP is the neutrahno, the hadrons produced in 
the neutrahno decays typically dissociate the primordial elements [4], yielding abun- 
dances in conflict with observations. On the other hand, when the NLSP is a charged 
particle, X~, the formation of the bound state ("^HcX^) catalyzes the production of 
^Li [5] leading to an overproduction of ^Li by a factor 300 — 600 [6]. 

Although the scenario depicted above is the most widely studied, it is not the only 
possibility. Indeed, several alternatives have been proposed that yield a thermal history 
of the Universe consistent with the observed relic density of dark matter, successful 
leptogenesis and successful Big Bang nucleosynthesis. For instance, in some specific 
supersymmetric models the NLSP can be a sneutrino [7] or a stop [8], whose late 
decays do not substantially affect the predictions of Big Bang nucleosynthesis. Another 
possibility is to assume some amount of entropy production after NLSP decoupling, 
which dilutes the NLSP abundance [9]. Finally, if i?-parity is not exactly conserved, the 
NLSP can decay into two Standard model particles well before the onset of Big Bang 
nucleosynthesis, avoiding the BBN constraints altogether [10]. This is the scenario that 
we will adopt in this paper. 



When i?-parity is not imposed, the superpotential of the Minimal Supersymmetric 
Standard Model (MSSM) reads [11] 



where Wr^ is the familiar superpotential with conserved i?-parity. Present laboratory 

experiments very severely restrict the size of the i?-parity breaking couplings. For in- 
stance, when the soft masses arc ~ 100 GeV, proton stability requires A'^]^^,A"]^^, ^ 10"^''', 
and the non-observation of the lepton flavor violating process fiTi eTi requires 
^ik2^'ku ~ 4 X 10~^, for /c = 1, 2, 3. An exhaustive list of the laboratory constraints on 
the it!-parity violating couplings can be found in [12]. 
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In addition to the laboratory upper bounds, there also exists an allowed window 
for the -R-parity violating Yukawa couplings stemming from cosmology. If the im- 
parity violating interactions had been in thermal equilibrium before the electroweak 
phase transition, any preexisting baryon or lepton asymmetry would have been erased. 
Therefore, successful leptogenesis can only be achieved if the out-of-equilibrium con- 
dition A, A', A" < 10~^ is satisfied [13]. These bounds are sufficient but not necessary 
conditions and could be relaxed for some specific flavor structures. On the other hand, 
successful Big Bang nucleosynthesis is guaranteed if the NLSP hfetime is shorter than 
~ 10^ s, which yields a lower bound on the i?-parity breaking Yukawa couplings. For 
instance, when the NLSP is a right-handed stau, it can decay via tr — > iii/r through 
the coupling A323 with lifetime 



Therefore, even a tiny amount of i?-parity violation, A323 ^ 10~ , is enough to deplete 
the population of stau NLSPs at the time of BEN down to harmless levels [10]. A 
similar argument applies for the case of a neutralino NLSP with analogous conclusions. 

When i?-parity is not exactly conserved, the gravitino LSP is no longer stable. 
Nevertheless, the gravitino decay rate is doubly suppressed by the Planck mass and by 
the small i?-parity violation, yielding [14,15] 



Therefore, for the range of it!-parity violating couplings favored by cosmology, 10~ < 
A, A' < 10~^, the gravitino lifetime ranges between 10^^ and 10^^ s for 777.3/2 = 100 GeV, 
which exceeds the age of the Universe by many orders of magnitude. Hence, even 
though the gravitino is not absolutely stable, it is stable enough to constitute a viable 
candidate for the dark matter of the Universe, while preserving the attractive features 
of the standard Big Bang nucleosynthesis scenario and thermal leptogenesis. 

Interestingly, the gravitino decay products could be observed as a contribution 
to the flux of cosmic gamma rays, positrons, antiprotons and neutrinos, opening the 
possibility of indirect detection of gravitino dark matter. We computed in [16] the 
gamma ray spectrum from gravitino decay, and we found that the anomaly in the 
extragalactic gamma ray flux reported by Strong et al. between 2-10 GeV [17] in the 
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EGRET observations [18] could be qualitatively explained by the decay of gravitino 
dark matter with a mass of 7/13/2 — 150 GcV and a lifetime of r3/2 — 1-3 x 10^^ s^. The 
expected anisotropy in the diffuse gamma ray flux was also found to be consistent with 
the EGRET observations. 

Motivated by this result, in this paper we compute the predicted fluxes of positrons 
and antiprotons from gravitino decay for the same set of parameters, as independent 
tests of this scenario^. The flux of positrons has been measured by a series of experi- 
ments: HEAT [21], CAPRICE [22], MASS [23] and AMS-01 [24]. Clearly, if gravitino 
decay is the explanation for the extragalactic EGRET anomaly, our predicted positron 
flux should not exceed the measured one. Although the measurements still suffer from 
large uncertainties, it is intriguing that they seem to point to an excess of positrons 
at energies larger than 7 GeV, which is precisely the energy range where we expect 
a contribution of positrons from gravitino decay. On the other hand, the antiproton 
flux has been measured by BESS [25], IMAX [26] and WiZard/ CAPRICE [27]. The 
measurements do not show any deviation from the predictions by conventional astro- 
physical models of spallation of cosmic rays on the Milky Way disk. Therefore, the 
viability of our scenario requires that the total antiproton flux lie below the astrophysi- 
cal background. Future antimatter experiments such as PAMELA [28] or AMS-02 [29] 
will provide very precise measurements of the spectra of positrons and antiprotons and 
will provide important constraints on the scenario of decaying dark matter. 

In Section 2 we derive the source term for positrons and antiprotons from the decay 
of gravitinos in the Milky Way halo. Both species, being electrically charged, propagate 
through the halo in a complicated way that we simulate by means of a conventional 
diffusion model. In section 3 we discuss and solve the diffusion equation for positrons 
and antiprotons. Finally, we present our conclusions and an outlook in section 4. 

2 Source Term 

We will assume that the Milky Way dark matter halo is populated by gravitinos with 
mass 777.3/2, their distribution following a density profile p(r), where r denotes the 

^The lifetime quoted in [16], T3/2 — 1-6 x 10^^ s, corresponds to a local halo density pq = 
0.38 GcV/cm^. Here we have adopted the more standard value = O.SOGcV/cm^ [19]. Note 
that the flux of particles from gravitino decay is proportional to /Co/t3/2; therefore the uncertainty in 
the value of the local halo density translates into an uncertainty in the gravitino lifetime. 

^The predictions for the neutrino flux will be presented elsewhere [20]. 
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Halo model 


a /3 7 Tc (kpc) 


Navarro, Frenk, White [30] 
Isothermal 
Moore [31] 


13 1 20 
2 2 3.5 
1.5 3 1.5 28 



Table 1: Parameters characterizing some commonly used halo models. 



position with respect to the center of the Galaxy. The dark matter distribution is 

usually parametrized as a spherically symmetric profile 

^('^) = (r/re)^[l + (r/re)"](^-^)/" ' 

where r = |r| and the parameters a, /3, 7 and Tc are listed in Table 1 for some commonly 
used halo models. Finally, po is a parameter that is adjusted to yield a local halo density 
of p{rQ) ~ 0.30GeV/cm^ [19], with ~ 8.5 kpc being the distance of the Sun to the 
Galactic center. 

Gravitinos at r eventually decay with lifetime producing antimatter at a rate 
per unit energy and unit volume given by 

Q^E,r, = ^^§, (7) 

where dN/ dE is the energy spectrum of antiparticles produced in the decay. 

If the gravitino is lighter than the gauge bosons, the main decay channel is 
"^3/2 which does not produce antimatter. On the other hand, if it is heavier than 

the gauge bosons, the decay channels iIjs/2 W^i"^ and 1IJ3/2 ai'c kincmatically 

accessible and produce antimatter"'. Namely, the process ipz/2 W~e^ produces a 
high-energy monoenergetic positron. On the other hand, the antimuon and antitau 
produced in the processes ■^3/2 — > W" p& ,W~t'^ generate a continuous spectrum of 
positrons in their decays. Lastly, the main contribution comes from the fragmenta- 
tion of the and the gauge bosons, which produce a continuous spectrum of 
positrons (mainly from 7r+ decay) and antiprotons that we have obtained using the 
event generator PYTHIA 6.4 [32]. Thus, the total energy spectrum of antiparticles 



^The size of the i?-parity violating couplings for the third generation is expected to be larger than 
for the first and second generations. Consequently, the charged lepton and the neutrino produced in 
these decays are expected to have predominantly tau flavor. 



5 



reads: 

dN dN^^ d N^" 

^ C BR(V.3/2 - ^ + BR(V^3/2 - Z%) ^ . (8) 

where dN^^/dE and dN^'^/dE denote the energy spectra of antiparticles produced 
in the indicated decay channel, which depend only on the gravitino mass through the 
total available energy. 

On the other hand, the branching ratios of the relevant decay channels can be 
straightforwardly computed from the expressions for the gravitino decay rates in [16]. 
The result is: 



BR(^3/2 ^ Wi) 



0\TT- |2 f( Mw - 



Mz 



\U- |2 f(Mz_) 



where f{x) = 1 — + ^x^ and U;yu, U-^^, U^y denote the mixings photino-neutrino, 
charged wino-charged lepton and zino-neutrino, respectively, which satisfy the following 
relations: 

(M2 — Ml) sin 9y/ cos 6w 



\U~ ~ 



Ml cos2 Ow + M2 sin^ Ow 



\U~zv\ > 



\Uwe\ - v2cos% — \U^J . (10) 

In this expression. Mi and M2 are the U{1)y and SU{2)i gaugino masses at low 
energies and 9\y is the weak mixing angle. It is commonly assumed that the gaugino 
masses unify at the Grand Unified Scale, Mx = 2 x 10^*^ GcV. Under this assumption, 
the ratio between the gaugino masses at low energies is predicted to be M2/M1 ~ 1.9, 
which yields 

: |[/^J : \U^^\ ~ 1 : 3.2 : 3.5 . (11) 

Therefore, the gravitino branching ratios in the different decay modes depend only on 
the gravitino mass (see Table 2). 

We conclude that, under the assumption of gaugino mass universality, the total 

energy spectrum of antiparticles from gravitino decay, dN/dE, depends exclusively on 
the gravitino mass. This makes our scenario very predictive: for a given halo model the 
source term Q{E,r) depends on only two unknown parameters, namely the gravitino 
mass and the gravitino lifetime; the former determines the spectral shape of the source 
function and the latter the normalization. 
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^3/2 


BR('03/2 iv) 


BR(V3/2 ^ Wl) 


BR(V3/2 ^ Z^v) 


10 GeV 


1 








85 GeV 


0.66 


0.34 





100 GeV 


0.16 


0.76 


0.08 


150 GcV 


0.05 


0.71 


0.24 


250 GeV 


0.03 


0.69 


0.28 



Table 2: Branching ratios for gravitino decay in different i?-parity violating channels for 
different gravitino masses. 

3 Antimatter Propagation in the Galaxy 

Antimatter propagation in the Milky Way is commonly described by a stationary two- 
zone diffusion model with cylindrical boundary conditions [33]. Under this approxima- 
tion, the number density of antiparticlcs per unit kinetic energy, f{T,f,t), satisfies the 
following transport equation, which applies both for positrons and antiprotons: 

= ^ = V • [K{T, f) V/] + ^[b{T, r)f] - V ■ [V,{r)f] - 2h6{z)T,^J + Q{T, r) . (12) 

The boundary conditions require the solution f{T, f, t) to vanish at the boundary of 
the diffusion zone, which is approximated by a cylinder with half-height L — 1 — 15 kpc 
and radius R — 20 kpc. 

The first term on the right-hand side of the transport equation is the diffusion term, 
which accounts for the propagation through the tangled Galactic magnetic field. The 
diffusion coefficient K(T, r) is assumed to be constant throughout the diffusion zone 
and is parametrized by: 

K{T) = Kopn^ (13) 

where (5 — v/c and TZ is the rigidity of the particle, which is defined as the momentum 
in GeV per unit charge, TZ = p{GeV)/Z. The normalization Kq and the spectral index 
5 of the diffusion coefficient are related to the properties of the interstellar medium 
and can be determined from the fiux measurements of other cosmic ray species, mainly 
from the Boron to Carbon (B/C) ratio [34]. The second term accounts for energy losses 
due to inverse Compton scattering on starlight or the cosmic microwave background, 
synchrotron radiation and ionization. The third term is the convection term, which 
accounts for the drift of charged particles away from the disk induced by the Milky 
Way's Galactic wind. It has axial direction and is also assumed to be constant inside 
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the diffusion region: Vc{r) = Vc sign(z) k. The fourth term accounts for antimatter 
annihilation with rate rann, when it interacts with ordinary matter in the Galactic 
disk, which is assumed to be an infinitely thin disk with half-width h = 100 pc. Lastly, 
Q{T,r) is the source term of positrons or antiprotons which was derived in Section 2. 
In this equation, reacceleration efi^ects and non-annihilating interactions of antimatter 
in the Galactic disk have been neglected. 

The solution of the transport equation at the Solar System, r — Tq, z — 0, can be 
formally expressed by the convolution 

1 /-Tmax dN(T'\ 

f{T) = -—— / dT'G{T, T')-^ , (14) 

where Tmgx = ^77,3/2 for the case of the positrons and Tmax = "^3/2 — trip for the an- 
tiprotons. The solution is thus factorized into two parts. The first part, given by the 
Green's function G{T,T'), encodes all of the information about the astrophysics (such 
as the details of the halo profile and the complicated propagation of antiparticles in the 
Galaxy) and is universal for any decaying dark matter candidate. The remaining part 
depends exclusively on the nature and properties of the decaying dark matter candi- 
date, namely the mass, the lifetime and the energy spectrum of antiparticles produced 
in the decay. 

Finally, the flux of primary antiparticles at the Solar System from dark matter 
decay is given by: 

^pHm(y) = ^fiT), (15) 

where v is the velocity of the antimatter particle. 

In the scenario we are considering the gravitino mass and lifetime are constrained by 
requiring a qualitatively good agreement of the predicted extragalactic gamma ray flux 
with the EGRET data: m3/2 = 150 GeV and T3/2 = 1-3 x 10^^ s [16]. On the other hand, 
the energy spectrum of antiparticles, dN/dT, is determined by the well- understood 
physics of fragmentation. Therefore, the only uncertainties in the computation of the 
antimatter fluxes stem from the determination of the Green's function, i.e. from the 
uncertainties in the propagation parameters and the halo profile. As we will see, the 
uncertainties in the precise shape of the halo profile are not crucial for the determination 
of the primary antimatter fiuxes, since the Earth receives only antimatter created 
within a few kpc from the Sun, where the different halo profiles are very similar. 
On the other hand, the uncertainties in the propagation parameters can substantially 
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change the predictions for the antimatter fluxes, even by two orders of magnitude for 

the antiproton flux. 

The reason for this large uncertainty is a correlation among the diffusion parameters 
and the size of the diffusion zone. Secondary cosmic rays are produced by spallation 
of primary cosmic rays in the Galactic disk. Therefore, the measurement of primary 
and secondary cosmic ray fluxes (particularly the Boron to Carbon ratio) provides 
information about the diffusive properties of the interstellar medium. Unfortunately, 
there exist degeneracies in the determination of the diffusion parameters. For instance, 
an increase in the size of the diffusion zone, which allows for a longer propagation 
time of cosmic rays inside the diffusion zone before escaping, can be compensated by a 
simultaneous increase of the diffusion coefficient, which facilitates a faster diffusion of 
cosmic rays away from the Galactic disk. However, this degeneracy does not hold for 
the antimatter fluxes from decaying dark matter, since antimatter is not only produced 
in the Galactic disk, but in the whole dark matter halo. Therefore, an increase in the 
size of the diffusion zone translates into an increase in the number of injected primary 
antiparticles, which is not compensated by the simultaneous increase of the diffusion 
coefficient. As a result, the antimatter fluxes from decaying dark matter can vary 
substantially for the range of astrophysical parameters which successfully reproduce 
the secondary cosmic ray fluxes. The ranges of the astrophysical parameters that are 
consistent with the B/C ratio and that produce the maximal, median and minimal 
positron and antiproton fluxes are listed in Tables 3 and 5 [35]. 

Positrons and antiprotons have different properties and their respective transport 
equations can be approximated by different limits of Eq. (12), thus allowing simple 
analytic solutions. Let us discuss each case separately. 

3.1 Positron Flux 

For the case of the positrons. Galactic convection and annihilations in the disk can be 
neglected in the transport equation, which is then simplifled to: 

V • [K(T, r) V/e+] + ^[b{T, r)/e+] + Q(T, r) = , (16) 

where the rate of energy loss, b{T,r), is assumed to be a spatially constant function 
parametrized by b{T) = with Tq = 1 GeV and te = 10^^ s. 
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Model 


6 


Ko (kpcVMyr) 


L (kpc) 


M2 


0.55 


0.00595 


1 


MED 


0.70 


0.0112 


4 


Ml 


0.46 


0.0765 


15 



Table 3: Astrophysical parameters compatible with the B/C ratio that yield the minimum 
(M2), median (MED) and maximal (Ml) flux of positrons. 



The solution to this equation is formally given by the convolution Eq. (14). The 
explicit form of the Green's function is [36] 

GeAT, T') = £ ^ BUT. T')Jo (Cn^) sin (^) . (17) 

where Jq is the zeroth-order Bessel function of the first kind, whose successive zeros 
are denoted by On the other hand, 



s„^(r,r') = ^C7„^exp 



with 

2 fR 




TJ VTn 



(18) 



2L ^ ' 



(19) 



where Ji is the first-order Bessel function. 

The Green's function can be well approximated by the following interpolating func- 
tion, which is valid for any decaying dark matter particle: 

1 ni6 

(T, T') ~ i^ea+6{T^-i-T'^-i)^( j^/ _ T) cm-3 s , (20) 

where T and T' are expressed in units of GeV. The coefficients a and b can be found in 
Table 4 for the NEW profile and the different diffusion models hsted in Table 3. This 
approximation works better than a 15-20% over the whole range of energies. We find 
numerically that the Green's function is not very sensitive to the choice of the halo 
profile; therefore the corresponding coefficients can be well approximated by Table 4. 

The interstellar positron flux from gravitino decay can be computed from Eqs. (14) 
and (15), the result being: 

$P^i-(T) = P'' dT'G,+ (T,T') '^^'tF'^ . (21) 

^ ' 47rm3/2r3/2 ^ ^ ' ' dT' ^ ' 

10 



model 


a h 


M2 
MED 
Ml 


-0.9716 -10.012 
-1.0203 -1.4493 
-0.9809 -1.1456 



Table 4: Coefficients of the interpolating function Eq. (20) for the positron Green's function, 
assuming a NFW halo profile and for the different diffusion models in Table 3. 



Secondary positrons 



Moore et al. ^,-'''^'r.r-' 

,>'^^r>*-*''"' NFW 
Isothermal 





Secondary positrons 






Ml 






MED 






M2 


\ ^ 
I, 



10 
T[QeV] 



10 
T[QeV] 



Figure 1: Interstellar positron flux from the decay of gravitinos with 771,3/2 — 150 GcV and 
■^3/2 — 1-3 X 10^^ s. In the left plot we assume the M2 diffusion model (see Table 3) and we 
study the sensitivity of the positron flux to various halo profiles. On the other hand, in the 
right plot we assume a NFW halo profile and we study the sensitivity of the positron flux to 
the diffusion model. We also show for comparison the secondary positron flux from spallation 
of cosmic rays on the Galactic disk. 

We show in Fig. 1 the predicted interstellar positron flux from gravitino decay for 
different halo profiles (left plot) and for different diffusion models (right plot). As 
expected, the dependence of the positron flux on the choice of the halo model is quite 
weak. On the other hand, the dependence on the diffusion model is important only at 
low energies, where the signal lies well below the background. At energies where the 
contribution to the total positron flux from gravitino decay can be visible, T > 7 GeV, 
the choice of the diffusion model only changes the primary positron flux by a factor 
2-3. 

Rather than measuring the positron flux, most experiments measure the positron 
fraction, $e+/($e- + ^e+)) since most sources of systematic error, such as detector 
acceptance or trigger efficiency, cancel out when computing the ratio of particle fluxes. 
Furthermore, the effects of solar modulation, which are important in computing the 
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positron flux at the top of the atmosphere below 10 GeV, also cancel out in the positron 
fraction when solar modulation is assumed to be charge-sign independent. In addition 
to the primary positron flux from gravitino decay there exists a secondary positron 
flux originating from the collision of primary protons and other nuclei on the inter- 
stellar medium, which constitutes the background to our signal. For the background 
fluxes of primary and secondary electrons, as well as secondary positrons, we use the 
parametrizations obtained in [37] from detailed computer simulations of cosmic ray 
propagation [38]: 

^r(^) - i + iiT'^+lar^-i^ (GeV-WV^sr-^) , (22) 

*-(^) = l+65Qr'3 + 16QQr^^^ (GeV-cm-Vs.-') . (24) 
where T is expressed in units of GeV. Then, the positron fraction reads: 

PF(T) = : l^J + '^e+l^J /25) 

where following [37,39] we have left the normalization of the primary electron flux as 
a free parameter, k, to be fltted in order to match the observations of the positron 
fraction. When there is no primary source of positrons, the positron fraction is best 
fltted for k = 0.88 [39]. 

We show in Fig. 2 the positron fraction for different halo proflles (left plot) and for 
different diffusion models (right plot). In accordance with the results for the primary 
positron flux, the dependence of the positron fraction on the halo model is very weak. 
Furthermore, the mild dependence of the primary positron flux on the choice of the 
diffusion model becomes even milder when computing the positron fraction. The reason 
for this is double: flrstly, the primary positron flux is never much larger than the 
secondary positron flux, and secondly, the dependence on the choice of diffusion model 
is partially absorbed by the normalization of the primary electron flux that we have 
left as a free parameter. Note also that the M2 model, which produces the minimal 
primary positron flux, yields the most prominent bump in the positron fraction. This 
"inversion" is again a consequence of having left the normalization of the primary 
electron flux as a free parameter. In order to reproduce the measured positron fraction 
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Figure 2: Same as Fig. 1, but for the positron fraction. 

at low energies, the normahzation of the primary electron flux k has to be smaller in 
the M2 model than in the MED and Ml (the precise values are k = 1.07, 1.28, 1.29 
for the M2, MED and Ml model respectively). Then, with the primary positron flux 
being comparable for all the diffusion models at energies above ~ 10 GeV, the smaller 
value of k for the M2 model yields a larger positron flux in this energy range than for 
the Ml and MED models. 

In conclusion, we find that gravitino parameters which predict a departure from a 
simple power law in the extragalactic gamma ray spectrum at energies above 2 GeV (as 
observed by EGRET), inevitably predict a bump in the positron fraction at energies 
above 7 GeV (as observed by HEAT). Furthermore, the presence of this feature is 
not very sensitive to the many astrophysical uncertainties. This remarkable result 
holds not only for the scenario of gravitino dark matter with broken i?-parity, but also 
for any other scenario of decaying dark matter with lifetime ~ 10^^ s which decays 
predominantly into and/or gauge bosons with momentum ~ 50 GeV. 

3.2 Antiproton Flux 

The general transport equation, Eq. (12), can be simplified by taking into account that 
energy losses are negligible for antiprotons. Therefore, the transport equation for the 
antiproton density, fp(T,r,t), is simplified to: 

df- 

= l^f = ^ ■ (^(^' ^1^/?) - ^ ■ - 2M(z)r,,„/p- + Q(T, f) , (26) 
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Model 


6 


Ko (kpcVMyr) 


L (kpc) 


V;(km/s) 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



Table 5: Astrophysical parameters compatible with the B/C ratio that yield the minimal 
(MIN), median (MED) and maximal (MAX) flux of antiprotons. 



where the annihilation rate, Fann, is given by 

Tann = (^H + ^'/^nne)<J^''Vp . (27) 

In this expression it has been assumed that the annihilation cross section between an 
antiproton and a helium nucleus is related to the annihilation cross section between 
an antiproton and a proton by the simple geometrical factor 4^/^. On the other hand, 
~ 1 cm~^ is the number density of Hydrogen nuclei in the Milky Way disk, ~ 
0.07 riH the number density of Helium nuclei and app^ is the annihilation cross section, 
which is parametrized by [40]: 

aiin/^N _ / 661 (1 + 0.0115 T-°-^^^ - 0.948 T^-^^^^) mbarn , T < 15.5 GeV , , . 
^PP ~ I 36 mbarn , T > 15.5 GeV , ^ ^ 

Analogously to the positron case, the solution to the transport equation can be 
expressed as a convolution of the form Eq. (14). The analytic expression for the 
Green's function reads [41]: 

G,(T,r') ^ge.p(- J^) ^^^-^-l^^^^cr-l)s(T-r) , (29) 

where 

(30) 



and 



A{T) = 2hr^^{T) + V, + kSi{T)coth^^^, (31) 



S^{T) = 



4C^ 

' +4- (32) 
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model 


X y z 


MIN 
MED 
MAX 


-0.0537 0.7052 -0.1840 
1.8002 0.4099 -0.1343 
3.3602 -0.1438 -0.0403 



Table 6: Coefficients of tlie interpolating function Eq. (33) for the antiproton Green's function 
for the NFW halo profile. 



We find that the Green's function can be numerically approximated by the following 
interpolation function: 

Gp{T, T') ~ 10^^ ^x+ylnT+zln^T^^j., _ j.^ ^^-3 ^ ^ (33) 

which, again, is valid for any decaying dark matter particle. The coefficients x, y and z 
for the NFW profile can be found in Table 6 for the various diffusion models in Table 5. 
In this case the approximation is accurate to a 5-10%. As in the case of the positrons, 
the dependence of the Green's function on the halo model is fairly weak. 
The interstellar antiproton flux is then given by 



^f{T) ^ ^ ^ \ dT'Gp{T,T') ^ . (34) 

f ^ ^ 47rm3/2r3/2 Jo p\ ^ J ^j., \ ) 

However, this is not the antiproton flux measured by balloon or satellite experiments, 
which is affected by solar modulation. In the force field approximation [42] the effect of 
solar modulation can be included by applying the following simple formula that relates 
the antiproton flux at the top of the Earth's atmosphere and the interstellar antiproton 
flux [43]: 

^^™^^-|, 2m,ris + Ti J^p-^^^s)' (35) 
where Tis = Ttoa + with Tjs and Ttoa being the antiproton kinetic energies at 
the heliospheric boundary and at the top of the Earth's atmosphere, respectively, and 

being the solar modulation parameter, which varies between 500 MV and 1.3 GV 
over the eleven-year solar cycle. Since experiments are usually undertaken near solar 
minimum activity, we will choose (pp = 500 MV for our numerical analysis in order to 
compare our predicted flux with the collected data. 

We show in Fig. 3, left plot, the predicted antiproton flux for different halo models. 
As in the case of the positrons, the sensitivity of the primary antiproton flux to the 
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Figure 3: Same as Fig. 1, but for the primary antiproton flux at tlie top of the atmosphere. 
In the left plot the MIN diffusion model was assumed (see Table 5). 

choice of halo model is fairly mild. We also show in the right plot the predicted 
antiproton flux from gravitino decay for the diffusion models listed in Table 5. From 
the plot, the extreme sensitivity of the primary antiproton flux to the choice of the 
diffusion model is apparent: parameters that successfully reproduce the observed B/C 
ratio lead to antiproton fluxes that span over two orders of magnitude. For a wide range 
of propagation parameters, the total antiproton flux is well above the observations and 
thus our scenario is most likely excluded, in spite of all the simplifying assumptions in 
the diffusion model. However, the MIN model yields a primary flux that is below the 
measured flux and thus might be compatible with observations. 

We have analyzed more carefully the predictions for the MIN model by computing 
the total antiproton flux. The result is shown in Fig 4, where, for consistency, we 
have adopted as background the secondary antiproton flux calculated in [41] for the 
same MIN model. Although the primary antiproton flux is smaller than the measured 
one, the total antiproton flux is a factor of two above the observations. Nevertheless, 
in view of all the uncertainties that enter in the calculation of the antiproton flux, it 
might be premature to conclusively rule out the scenario of decaying gravitino dark 
matter. Namely, in addition to the uncertainties stemming from degeneracies in the 
diffusion parameters, there are also uncertainties from the nuclear cross sections and, to 
a lesser extent, uncertainties from the description of the interstellar medium and solar 
modulation (for a discussion of the various uncertainties see [41]). Furthermore, we used 
a simplifled diffusion model that neglects the effects of reacceleration, energy losses and 
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Figure 4: Contributions to the total antiproton flux in the MIN diffusion model. 

tertiary contributions. Therefore, there could be certain choices of parameters or more 
refined diffusion models where the total antiproton flux is consistent with experiments^. 

4 Conclusions and Outlook 

In this paper we have calculated the positron and antiproton fluxes from gravitino dark 

matter decay. The source term merely depends on two parameters, the gravitino mass 

and the gravitino lifetime, rendering a very predictive scenario from the particle physics 

point of view. The main uncertainties arise from the astrophysics, namely from our 

ignorance of the precise shape of the halo profile and especially from the degeneracies 

in the determination of the diffusion parameters. 

By requiring a qualitatively good agreement of the predicted extragalactic gamma 

ray flux to the EGRET data, we have fixed 7713/2 = 150 GeV and r3/2 = 1.3 x 10^^ s. 

This choice of parameters completely fixes the source term, and the only indeterminacy 

in the computation of the antimatter fluxes stems from the unknown astrophysical 

parameters. Remarkably, with independence of the astrophysical uncertainties, we 

predict a bump in the positron fraction at energies above 7 GeV, in agreement with 

the HEAT observations. On the other hand, the predicted antiproton flux tends to be 

too large, although for certain choices of the propagation parameters the predicted flux 

might also be in agreement with observations. 

^Some works have reported a deficit in the predicted secondary antiprotons compared to the obser- 
vations and argued that this deficit could be connected with a contribution of primary antiprotons [44] . 
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The main conclusion of this paper is summarized in the three plots in Fig. 5. There, 
we show the predicted extragalactic gamma ray flux, positron fraction and antiproton 
flux compared to the EGRET, HEAT and BESS data respectively, for 171^/2 = 150 GeV, 
T3/2 = 1.3 X 10^^ s and the MIN diffusion model in Table 5. It is intriguing that for 
this diffusion model the scenario of gravitino dark matter with broken i?-parity can 
qualitatively explain the anomalies observed in the extragalactic gamma ray flux and 
the positron fraction in a very natural way. It should also be stressed that this scenario 
was not devised to explain the anomalies in the cosmic ray fluxes, but to reconcile the 
clashing paradigms of supersymmetric dark matter, thermal leptogenesis and Big Bang 
nucleosynthesis. At the same time, we flnd that the total antiproton flux is slightly 
larger than the observed one. However, given all the uncertainties that enter in the 
calculation of the antiproton flux, it might be premature to rule out the present scenario 
on the basis of this small excess. 

To conclude on the phenomenological viability of this scenario, it would be worth- 
while to elaborate on the propagation of positrons and especially antiprotons from 
gravitino decay by going beyond the simplifled diffusion model used in this paper [45] . 
On the experimental side, the upcoming gamma ray experiment GLAST and the an- 
timatter experiment PAMELA will provide in the near future measurements of the 
cosmic ray fluxes with unprecedented accuracy, thus providing invaluable information 
about the scenario of decaying gravitino dark matter. 

Finally, we would like to mention that these results are not peculiar to the scenario of 
gravitino dark matter with broken i?-parity. The characteristic feature of our scenario 
is that the dark matter decays at late times into gauge bosons, eventually producing 
photons, positrons and antiprotons in the fragmentation. Therefore, similar signatures 
can be expected for other decaying dark matter particles that couple to the and/or 
gauge bosons. 
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Figure 5: Summary of the signatures of gravitino dark matter decay in the extragalactic 
gamma ray flux (top), the positron fraction (bottom left) and the antiproton flux (bottom 
right), compared to the EGRET, HEAT and BESS data respectively. In these plots, we have 
adopted the MIN diffusion model (see Table 5), m3/2 — 150 GeV and r3/2 — 1.3 x 10^^ s. 
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